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Abstract 

We study the probability of formation of ferromagnetic string in the antiferromag- 
netic spin-1/2 XXZ chain. We show that in the limit of long strings with weak 
magnetization per site the bosonization technique can be used to address the prob- 
lem. At zero temperature the obtained probability is Gaussian as a function of the 
length of the string. At finite but low temperature there is a crossover from the 
Gaussian behavior at intermediate lengths of strings to the exponential decay for 
very long strings. Although the weak magnetization per site is a necessary small 
parameter justifying our results, the extrapolation of obtained results to the case of 
maximally ferromagnetic strings is in qualitative agreement with known numerics 
and exact results. The effect of an external magnetic field on the probability of 
formation of ferromagnetic strings is also studied. 

Key words: integrable models, spin chains, correlation functions, bosonization, 
emptiness formation probability, asymptotic behavior 



1 Introduction 



Antiferromagnetic spin-1/2 XXZ chain is the first integrable model known to 
be solvable by means of Bethe Ansatz [1,2]. It is described by the Hamiltonian 
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where J > is an exchange integral, A is an exchange anisotropy, and sum- 
mation is taken over lattice sites of an infinite spin chain. In our notations the 
particular values A = +1,0, —1 correspond to isotropic antiferromagnet, XY 
magnet, and isotropic ferromagnet respectively. It is well known [3-5] that for 
— 1 < A < 1 (so-called critical regime) the spectrum of excitations around 
the ground state is gapless and spin-spin correlation functions decay with the 
distance as a power law. In this paper we consider only the critical regime 
-1 < A < 1. 

Although thermodynamic properties of such a spin chain are well understood 
[4], the exact calculation of correlation functions is a long-standing hard prob- 
lem (see [3] and references therein). Correlation functions in the model (1) (as 
well as in other integrable models) in thermodynamic limit can be represented 
as combinations of minors and determinant of Fredholm operator of a special 
form [integrable integral operator]. From the point of view of this determinant 
representation the simplest and the most fundamental object is not a correla- 
tion function of local fields which is given by minors of Fredholm operator, but 
rather a probability P n of spontaneous formation of magnetization on some 
space interval (ferromagnetic string) of the length n which is given by the 
determinant of that operator. This probability was first introduced in [6] (see 
also [3]) as "emptiness formation probability" or "probability of formation of 
ferromagnetic string" (PFFS). It is given by 



where (3 = 1/T is an inverse temperature and Z — Tr je is a partition 
function of a spin chain (1). 

Another approach to correlation functions was developed in Kyoto[9]. In this 
approach correlations were represented as multiple integrals. Analysis of these 
integrals also shows that PFFS is the simplest correlation function. 

Unfortunately, both determinant and multiple integral representations of cor- 
relation functions grow extremely complex with the increase of the length of 
a ferromagnetic string. So far the PFFS was calculated exactly by means of 
multiple integral representation only up to the length n = 5[10] . Numeric 
results known for this probability suggest [10] that the asymptotic behavior 
is Gaussian P n ~ e~ an at zero temperature. At the moment there are only 
two exact results known about the large n behavior of P n . In both cases the 
Gaussian behavior is obtained for P n at large n. The first case is an XY spin 
chain (A = 0) which is equivalent to free lattice fermions. In this case an 
alternative determinant representation for P n is available which allows to cal- 
culate PFFS exactly for any length of the string[ll]. The large n asymptote of 
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PFFS is Gaussian (see [11]) P n ~ n~ 1 / 4 (y/2)~ n2 . Another important example 
is XXZ spin chain at A = 1/2 . In this case PFFS is also known explicitly [12] 
(although it is still a conjecture). In this case the large n behavior is again 
Gaussian P n ~ n" 5 / 36 ( v / 3/2) 3n2 . 

On the other hand there is a relatively simple and very powerful method of 
calculation of asymptote of correlation functions for XXZ spin chains. It is 
a combination of bosonization technique and Bethe Ansatz (for review and 
references see [5,7,8]). Namely, at large distances (in the continuum limit) 
the XXZ spin chain is equivalent to a free boson theory with parameters (e.g., 
compactification radius of bosonic field) which can be determined exactly from 
Bethe Ansatz results for thermodynamics of the chain. This approach was 
successfully used to calculate the long distance asymptotes of local spin-spin 
correlation functions in spin chains [5,7,8]. 

In this paper we apply the bosonization technique to study PFFS on the 
background of an antiferromagnetic spin chain in the limit of the large length 
n of the string. 

First, we give a qualitative argument in favor of the Gaussian decay of PFFS at 
zero temperature. We are looking for optimal space-time configuration of spins 
(instanton) which costs the least action but gives the maximal magnetization 
of the string of the length n at some moment of time. We use two assumptions. 
First, we assume that the effective theory describing spin chain is Lorentz 
invariant at large distances (it is actually anisotropic but scaling of space and 
time is the same). Second, the extra cost in action of having ferromagnetic 
state is constant per unit space-time area. This is a consequence of the first 
assumption and the fact that the cost of energy is constant given by exchange 
integral per lattice spacing. We conclude from these two assumptions that the 
optimal instanton has a shape of a circular (elliptic) droplet of the area ~ n 2 
with constant action cost per unit area. This immediately gives the minimal 
action ~ n 2 and P n ~ e~ an at zero temperature. The fluctuations around the 
found optimal configuration can only change pre-exponential factor but not the 
value of exponent itself. At finite temperature we expect similar behavior when 
the length of the string is much smaller than an inverse temperature (properly 
rescaled to spatial units). On the other hand, in the limit of long strings one 
of n's in the action of instanton is replaced by the inverse temperature (3 
( "time" is finite and instanton has a size n x (3. This will produce P n ~ e~ in 
- well known result [3] with 7 proportional to the density of free energy of an 
antiferromagnetic state. 

Unfortunately, bosonization is not directly applicable to finding PFFS — the 
probability of formation of maximally ferromagnetic string. The reason is that 
the maximal magnetization corresponds to the large (compared to an inverse 
lattice spacing) gradient of the bosonic field which violates the continuum ap- 
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proximation. It turns out, however, that one can use this technique to find 
the probability of formation of weakly ferromagnetic strings. For weakly fer- 
romagnetic strings our calculations justify the qualitative argument in favor 
of Gaussian decay given above and the obtained results are in qualitative 
agreement with numerics [10,13] and exact results [11,12] for the strings with 
maximal magnetization. 

In Sec. 2 we briefly review the well known results of bosonization of XXZ spin 
chain. In Sec. 3 we present our calculation of the probability of formation of 
weakly ferromagnetic string. We summarize the results of this calculation in 
Sec. 4 and compare them with available numeric and exact results in Sec. 5. 
We relegated the calculation in the presence of an external magnetic field and 
the calculation for strings with uniform magnetization to Appendices A, B. 



2 Bosonization 



Bosonization of the continuum limit of XXZ chain gives the action (for review 
and references see [5,7]) 



S 



dxjdr 
o 



^-{d^f + V -{d x ct>f 



(3) 



where x = j is a continuous coordinate along the chain, and lattice spacing 
is taken as a unity a = 1, r is an imaginary time (h = 1), v is a spin wave 
velocity, and R is a radius of compactification of free Bose field <fi: <fi = <f)+2-n-R. 
The parameters R and v of the model (3) are known exactly from the Bethe 
ansatz solution of the model (1) [3] 



2irR 2 = 1 - - cos" 1 A, (4) 



7T 



cos 1 A 

In the following we putQ v — 1. 

S= J dx j dr-(d^) 2 . (6) 

-oo 



1 One can always restore the dependence on spin wave velocity replacing at the end 
of calculation f3 — > (3 v. 
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The original spin operators are given in terms of bosonic field as 



a] = -Lg + const- (-iy cos | 



const ■ cos ^- + C(— 1 Y 
R 



where 



(f)(r,x) = (j) L (z) + <j> R (z) 
0( r ,x) = (f) L (z) -<f> R {z) 

and z = x + it, z = x — it are complex coordinates. 

Let us consider a long segment of the spin chain of the length n 1. The total 
magnetization of the segment at some time r is given by m n (r) = J2j °~] = 
^ / dxd x (() = ^ [4>{n, r) — 0(0, r)] . The maximal magnetization (ferromag- 
netic string) corresponds to J2j &j — n- We are looking for the probability 
of spontaneous formation of the string of spins of the length n and average 
magnetization per site c (c = 1 corresponds to the maximal magnetization 
a z = 1 per site). We require therefore, 

[A0] = <j>(n) - 0(0) = cnnR. (7) 



Here we denoted the change of the bosonic field as [A0] to avoid a confusion 
with the Laplacian of the same field A0. Notice that the typical value of the 
gradient of along the string is d x <j) ~ ci\R. 

Bosonization is an essentially continuum approach which is valid only in the 
limit when (9 X ~ (9 r irR. This requirement gives the condition c < 1 
with the original problem of finding the probability of maximally ferromag- 
netic string to be beyond the range of applicability of bosonization approach. 
In the following we assume cC 1 and look for the probability of weakly ferro- 
magnetic strings. The following remark is in place. In contrast to the problem 
of maximally ferromagnetic string the problem of finding the probability of 
a weakly ferromagnetic string is not yet completely defined. Namely, one can 
search for the probability specifying the total magnetization of the string only 
(7) or one can fix the magnetization profile along the string a z (x) = cnRf(x) 
with arbitrary real function f(x) satisfying /(0) = and f(n) = n. In the fol- 
lowing we adopt the first interpretation (obviously it is equivalent to solving 
the problem with given f(x) and then averaging over all f(x)). The problem 
in the second interpretation can also be straightforwardly solved by presented 
methods. We consider the particularly interesting case of uniform magnetiza- 
tion in Appendix B. To use the bosonization method one should also require 
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n^$> 1, i.e., long strings. Finally the free Bose field approximation can be used 
only at large space-time scales where irrelevant (and marginal irrelevant for 
XXX case) operators coming from a z a z term renormalize to zero. It turns 
out, however, that this happens at scales of the order of 1 and imposes no 
additional restriction. 



3 Calculation 

The probability to have a ferromagnetic string is given by 

P n = i J £>05([A0] - cnnR)e- s ^, (8) 



where S[(j)} is given by (6) and Z — J D(pe s ^. We represent delta function 

as 

P„ = | / y I D<t>e- s W- iX ^- c ™ R \ (9) 



-oo 



The obtained integral in 4> is Gaussian and is determined exactly by its saddle 
point which can be obtained from 

A0 O = i\ \S 2 (z-n) -S 2 (z)} (10) 



with periodic boundary conditions in time 

(j)(z + 0) = 4>{z) + 2nRw, (11) 



where w is any integer number (winding number). The expression 5 2 (z — n) — 
8{x — n)5(r) denotes two-dimensional Dirac delta function. In this paper we 
consider only the limit of an infinitely long spin chain. Then, only topologically 
trivial w = sector will contribute to the probability because the action for 
the sectors with w ^ is proportional to the length of the chain. Therefore, 
we require 

<f>(z + iP) = <f>(z). (12) 
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3.1 Zero temperature 



Let us first consider the case of zero temperature f3 — > oo. Then (12) is irrele- 
vant (instead we require the decay of d^cf) at infinity) and the solution of (10) 
is given by 



^o(z) = ^-Re In- — - 

Z7T Z 



(13) 



and being substituted into (9) gives the Green's function of free boson at zero 
temperature (e-^e*^ )) or 



+oo 



" J 2vr 



(14) 



where we introduced new notation 



Q 



CTlR 



> 1 



(15) 



and r is some short distance cutoff. This integral is in turn determined by its 
saddle point 



2tt Qln^ 



(16) 



and is given by 



21n^ 

r 



Q 2 ln-2- 

e r 



(17) 



Let us now find a short distance cutoff r . On one hand it must be at least big- 
ger than 1 - the lattice scale. On the other hand we can trust our bosonization 
approach only if <9 M <C nR. From (13) we have on saddle point configuration 



2tt \z 



n 



This gradient diverges near the ends of the string. At points close to z = 
(similar for points close to z — n) we have 



2ir 



(18) 
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where r = \z\ {pi r = \z — n\ near the other end of the string). We obtain r 
from 

|9„^o| r =ro = -7T-— = T-R- ( 19 ) 

2% r 



Using (16) we obtain 



-°=nRr = -hi-. 20 
27r Q r 



We solve the latter equation by iterations and obtain 

cn ( / 1 
21UC- 1 V + VMc 



In (21) and in the following we have actually replaced O ( ^~r~ ) by O • 
Equation (21) gives the short distance cutoff for long strings n ^> c -1 lnc -1 . 
If strings are not so long one should use the lattice spacing cutoff r ~ 1. We 
have 



cn 

r = max < 1, 



21HC" 1 

In the following we consider only sufficiently long strings n ^> c -1 lnc" 1 and 
use (21). We find therefore, that bosonization approach is valid everywhere 
except for the vicinity of the ends of the string which is defined by the radius 
To from (21) (r <C n). Substituting the found value of r into (17) we obtain 
(we keep only exponential dependence) 

P n ~ e~ an2 , when n > c _1 In c _1 (22) 
"=0^ (23) 

Let us now estimate the corrections to our results coming from the vicinities 
(of the radius r ) of the ends of the string where bosonization approach is not 
applicable. First of all, we find that the optimal configuration O is given by 
(13), (16), and (21) as 

^» = Ql^ Reln 7^' < 24 » 



We obtain 

[A0 O ] w (j) (n - r ) - (f) (r ) « 2— 
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which satisfies (as expected) the boundary conditions (7) with the accuracy 
1 / lnc -1 . We estimate an error coming from the vicinity to the end of the string 
assuming the maximal gradient d x (po ~ ttR in that vicinity 0o( r o) ~ 0o(O) ~ 
rodx4>o ~ TrRr ~ n/(Q In c -1 ) ~ [A0] with [A0] from (7). We obtain that 
in the phase difference calculations one can neglect the vicinities of the strings 
with our usual accuracy ~ ln ^_ t . Secondly, let us estimate the contribution 
to the action coming from those vicinities. Again we use <9 M 0o ~ kB. and 
find correction to the action r^d^o) 2 ~ (ttRtq) 2 ~ q2^ c -i ~ ar ^ i n c- 1 • 
We conclude that with the accuracy ln ^_ t one can neglect the regions where 
bosonization breaks down. 

Therefore, in the case of zero temperature the probability to have a weakly 
ferromagnetic string of the length n ^> c -1 lnc" 1 ^> 1 in the antiferromagnetic 
XXZ spin chain is Gaussian (22), (23) and the small parameter justifying our 
bosonization calculation is 1/lnc" 1 . 



3.2 Finite temperature 



In this section we consider small but finite temperature T such that (3 1 
(in full units f3v ^> 1 and one can still use continuum limit and bosonization 
technique. The only thing which is modified in this case compared to the case 
of zero temperature is that one should solve the same problem (9) on a cylinder 
< r < j3 with periodic boundary conditions (12). 

Similarly to the case or zero temperature we obtain the equation (10) for the 
saddle point configuration. With boundary conditions (12) its solution is given 
by 

Mz) = ^Reln r \ (25) 

2ix smh jjZ 

Notice that in the limit n /3 the main contribution to (9) comes from 
distances \z\ ~ n <^ (3 and one can expand hyperbolic sines in (25) obtaining 
(13). Thus, in the limit of low temperatures we restore the results obtained 
for zero temperature. For the probability of the string (compare to eq.(14)) 
we obtain 

f d\ -A£ ln si"MW/3) + - A 2n 

P n — — e 2n nr o/P « , (26) 
J 2n 



where we again assumed r <^ [3 in treating short distance cutoff. One can 
easily recognize in the first exponent the finite temperature correlator of free 
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Bose fields ( e - a ^W e a ^°)). Integral over A is determined by the saddle point 

i\o n 



2tt Q In ™M™//3) 



(27) 



and we find (similarly to eq.(17)) 



I ~To~ sinh(7rn//3) 

P„ = e ln . (28) 

" 2 } n smh(7ro/j3) V / 

7rr //3 



The gradient of the field 



d x (p - id T (p = ^ 



7T 71 

coth — (z — n) — coth —2; 



(29) 



diverges near the ends of the string and we again introduce the short distance 

2/3 COU1 /3 



cutoff ro as the scale such that d x ((>o(ro) = — 14 coth ?ro = 7nR. In the following 



we show that r <C /3. Using this inequality we obtain d x (j) (r ) = — 2^ — ^-R 



and find 



nRr ° = "2^- (30) 



Short distance cutoff r is defined by (30) and (27) as 



n 

gin" 



TlRr-Q = ^-tt 775V. (31) 

s-\ 1 sinh(7rn//3) v ' 



In the following we are interested in sufficiently long strings n ^> c _1 lnc _1 
and sufficiently low temperatures (3 ^> c _1 mc _1 . It is easy to show that in 
these limits the solution of (31) satisfies r 1. Using (31) we rewrite (28) as 

Pn = ~^e • (32) 



Equations (32) and (31) give the probability of weak ferromagnetic string at 
finite temperature. 



We solve (31) by iterations and obtain 



ttRtq = -, —, — t^-n-. (33) 

Qln^Q^^) 
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Next iterations give corrections which are small by parameter 



In ttRQ 



sinh(7rn//3) 
im/P 



< 



In cT 1 



< 1. 



Using (33) we obtain the leading exponent from (32) 



n Q 2 ln(TTRQ- 

K ~ e v 



7vn/f3 



(34) 



Let us now check that the use of bosonization to derive (34) is justified. Firstly, 
for bosonization to be applicable the main contribution to the magnetization 
should come from the scales large than r . Namely, using (25), (27), and (31) 
we obtain 



[A0 O ] ss cf) (n - r ) - </> (r ) 



2n 



(35) 



and magnetization (change of 0o) coming from the scale larger than r and 
given by (35) is indeed the total magnetization of the string. On the other hand 
the contribution from segments [0, r ] and [n— r , n] to the total magnetization 
is of the order of nRr (it is determined by maximal gradient equal to r d x (f)o ~ 
7ri?r ) and we obtain from (35) the required consistency condition 



txRtq < 



n 

Q 



which together with already used assumption r (3 gives 



r min(cn, (3). 



(36) 



Secondly, calculating the optimal action (the value in the exponent of (32)) 
we neglected contributions from the r -vicinities of the ends of the string. We 
estimate that neglected action as r^d^) 2 ~ (nRr ) 2 with maximal gradient 
of the order of irR. It is easy to see that if this consistency condition is satisfied 
one can neglect this contribution as it is smaller than the action coming from 
larger distances (irRr ) 2 <^ n(rrRr )/Q. 

We conclude that the use of bosonization is legitimate if the consistency con- 
dition (36) is satisfied. It is easy to check that the solution (33) indeed satisfies 
this condition in the limit of sufficiently long n ^> c -1 lnc" 1 strings and suffi- 
ciently low /3 ^> c _1 lnc _1 temperature. Under these conditions the equation 
(34) gives the probability of spontaneous formation of a ferromagnetic string 
of an arbitrary length and at arbitrary temperature. Below we consider two 
important limiting cases. 
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3.2.1 Short strings: 1 < n < f lnc" 1 

Let us assume that ~-^7g^ nRQ. This is equivalent to taking not very 
long strings n <^ ^ lnc -1 . In this case we drop the combination — from 
(34) and (33) and obtain 



ro = 2^ (37) 



and 



2 



Equations (37-38) reproduce the results (21-23) we obtained in the limit of 
zero temperature. 



3.2.2 Very long strings: n ^> fine 1 

In the opposite limit of long strings ^~~r^ ^ nRQ (or n ^> f lnc -1 ) we 
replace In {11 RQ shl ^™' p) ) -> ^ an d obtain from (33,34) 

= £ (39) 



and 



_ 2/3n 

P n ~ e 'Q T . (40) 



In this limit the probability decays exponentially with the length of the string 
as it is expected [3]. 



4 Results 



We summarize the results of our calculation as 



H I^J g ^ 7T71//3 ) r*^ < 



e " a ™ 2 ' «=q^ n«flnc-\ 



(41) 
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where P n is a probability of finding the weakly ferromagnetic string of the 
length n and total magnetization m n = cn in an antiferromagnetic XXZ spin 
chain. The result is valid for small per site magnetization c< 1, long strings 
n 3> c _1 mc _1 3> 1 and at finite but sufficiently low temperatures 1/T = 
P 3> c^lncT 1 » 1 (f3 is measured in units of l/t> with v from (5)). The 
large parameter Q = 2/(cttR) ^> 1 with compactification radius R from (4). 
In the limit of zero temperature the actual small parameter justifying our 
calculations is l/lnc _1 <C 1. 

It is interesting to notice from (41) that the crossover from exponential to 
Gaussian decay of probability occurs not at the scale n ~ (3 but at much 
larger (for weakly ferromagnetic strings cC 1) scale n ~ ^lnc -1 . There is a 
simple reason for such a behavior. We fixed only the total magnetization of the 
string. It turns out that the magnetization is concentrated mostly at the ends 
of the string, namely, within the distance of nj lnc -1 to the ends of the string. 
Therefore, the total length of the string is irrelevant and the crossover takes 
place only when the new scale becomes of the order of (3, i.e., n ~ /3 In cT 1 . 

We also present in this section the results we obtain in Appendix B (see 
(B.9,B.ll) for the probability of formation of a weak ferromagnetic string 
with uniform magnetization 



The small parameter justifying our calculation in the limit of zero temperature 
is c 2 in this case. It is much smaller for cC 1 than the parameter 1/ lnc^ 1 for 
the case of free string — the string with only the total magnetization fixed. 

It is easy to generalize our results to the case when external magnetic field h is 
present. The results are still given by (41) or (42) with the replacement c — > c 
where c is given in (A. 7) (this replacement means that only the deviation from 
the background magnetization induced by magnetic field is relevant). Also, 
all other parameters v, R must be taken from the exact solution in non-zero 
magnetic field (see Appendix A for details). 



5 Discussion 

We studied the probability of formation of ferromagnetic string in XXZ spin 
chains. The method of continuum bosonization we used is not applicable to the 
case of a maximal magnetization of the string. However, we were able to obtain 
the behavior of this probability in the limit of a weakly ferromagnetic strings. 




(42) 
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In this limit we have a controllable approximation with a small parameter 
related to the small per site magnetization of the string. Our main results 
are presented in (41,42) and in Appendix A in the presence of an external 
magnetic field. 

These results: crossover from Gaussian behavior at low temperatures (short 
strings) to an exponential decay at high temperatures (long strings) are in 
qualitative agreement with numeric results [10,13] and known exact results 
[11,12]. We have also presented the qualitative argument in favor of such a 
behavior in Sec. 1. 

Let us now see how far off we are quantitatively. We will try to predict the 
rate of Gaussian decay a by (illegally) extending our results (42) obtained in 
the limit of a weakly ferromagnetic string c<l with a uniform magnetization 
to the case of maximally ferromagnetic string c = 1 (which obviously has a 
uniform magnetization). We choose (42) versus (41) for such an interpolation 
because we it to be more precise for not very small c because of a milder 
divergence of the optimal field configuration in the case of a string with a 
uniform magnetization (see Appendix B). 

We take c = 1 (or Q = 2/txR) in (42) and obtain a = n 3 R 2 /8 with the 
exact value of compactification radius R from (4). This result is qualitatively 
correct giving a = at R = which corresponds to the ferromagnetic limit 
A = — 1. In this limit we do expect a disappearance of a Gaussian decay of the 
probability of PFFS (the range A < — 1 corresponds to Ising ferromagnetic 
behavior). However, the exact numeric result is not justified^. 

We use a = n 3 R 2 /8 with (4) and write 

^e a = e^( 1 ^ cos " lA ), (43) 

2 

where new parameter 5 is introduced so that P n ~ 5~ n . In fact, one can expect 
any coefficient of the order of one instead of 7r 2 /16 in (43). In the following 
table we compare the rates of Gaussian decay of PFFS in XXZ spin chains 
calculated from (43) with two exact results known for A = 0, | [11,12]. 



2 

The rate of Gaussian decay P n ~ 5~ n . 


A 


5 (exact) 


5 from Eq.(43) 





1.4142... (= V2 [11]) 


1.36 


1/2 


1.5396... (= (^) 3 [12]) 


1.51 



2 We are grateful to S. Lukyanov for stressing this point to us. 
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We see from this table that the agreement is reasonably good even quantita- 
tively. The comparison of (43) with known numerical results[10,13] obtained 
for different values of anisotropy A also gives a good agreement (~ 10%). 

We conclude with the remark that in addition to giving an insight to the 
nature of a Gaussian decay of the probability of formation of ferromagnetic 
strings in spin models the problem we considered is interesting on its own. 
For example, there are known systems where the dopant electron propagation 
is disfavored by the background antiferromagnetic state. In such situations 
the electron propagation will be assisted by spontaneous formation of ferro- 
magnetic strings. The actual magnetization profile will be determined by the 
competition between the benefit for the electron propagation and the action 
cost for the necessary ferromagnetic fluctuation. These are the details of a 
specific model which will determine the magnetization profile of an optimal 
fluctuation and a weakly ferromagnetic string is one of the candidates. 
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[8] prior to publication. A.G.A. would like to thank the Alfred P. Sloan Foun- 
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A The effect of an external magnetic field 



In this Appendix we show how to generalize the results obtained in Sec. 3 in 
the presence of an external magnetic field. Magnetic field is introduced into 
XXZ model as 



N 



N 



(A.l) 



The bosonization of this model gives[5,7] 



+oo /3 

S = J dx J dr 

-oo o 



(A.2) 



where we again use units in which v — 1 (at the end of calculation h h/v). 
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We are going to minimize (A. 2) with the condition (7) on the total magneti- 
zation of the string. We notice that the change of variables 

= 0' + \x. (A.3) 

7lR 



transforms (A. 2) and (7) into (we omit primes and drop an additive constant) 



+oo ft 



S= J dx J dr l -{d^)\ (A.4) 

-oo 

9n 

[A0] = ^, (A.5) 



where 



Comparing (6,7) and (A.4-A.7) we notice that the only effect of an external 
magnetic field on the probability of ferromagnetic string (41) is just some 
effective renormalization of parameter c as well as the renormalization of R and 
v. The renormalization of the parameter c has a very simple physical reason. 
An external magnetic field results in the uniform background magnetization 
h/(7i 2 R 2 ). The effective renormalization (A. 7) means that for the probability 
of the string it is the deviation of magnetization from the equilibrium one 
which is relevant. The dependence of parameters v , R on a weakF] magnetic 
field is also known [3] from the Bethe Ansatz solution. 



B Ferromagnetic strings with uniform magnetization 



The probability of formation of a weak ferromagnetic string is given by the 
minimum of action (6) with corresponding boundary conditions on Bose field 
0(z). In Sec. 3 we used (7) as a boundary condition fixing, therefore, only the 
total magnetization of the string. It is straightforward to minimize (6) with 
an arbitrary given profile of magnetization along the string. To find such a 

3 The weakness of magnetic field h/ir 2 R 2 <C 1 (or in full units h/ir 2 R 2 v <C 1) is 
required so that the field transformation (A.3) is within the applicability of bosoniza- 
tion approach. 
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minimum one should solve a two-dimensional Laplace equation (obtained by 
variation of (6) with respect to 0) 



A0 = O (B.l) 

with periodic boundary conditions 

<P{z + i(3) = <P{z) (B.2) 

and magnetization profile defined by 

2 

(j>(x, r = 0) = —f(x), when < x < n. (B.3) 

Here as usual Q = 2/{ckR). f(x) is an arbitrary (real) function satisfying 
/(0) = 0, f(n) = n so that the total magnetization of the string is cn and 
magnetization profile of the string is given by a z (x) = cd x f(x). It is easy to 
check that the solution of Dirichlet problem (B.1-B.3) is given by 

2tt Ids f(s) /sinh^sinh^^\ * 

0( z ) = g- Re /- ^7 . £ ^ . (B.4) 

PQ I m sinh^^ \sinhf sinh / V ' 

Indeed, it obviously satisfies Laplace equation and periodic boundary condi- 
tions. For real z = x with < x < n the only real contribution to the integral 
can come from the half of the residue at s = x which gives (B.3). In the fol- 
lowing we restrict ourselves to the simplest case of a uniform magnetization 
f(x) = x along the string. In this case we have 

2x 

<f>(x, r = 0) = — , when < x < n. (B.5) 



In the limit of zero temperature (3 — > oo the solution of (B.l, B.5) is given by 
(one can either check it directly or obtain it from (B.4)) 

M*) = ^ Re ( z - ~ n )) • ( B - 6 ) 



We shift the position of the string for notational convenience to — n < x < n 
with n = n/2 and obtain 

fa(z) = 1 Re (z - v^ 2 - n 2 )) • (B.7) 
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One can recognize in (B.7) the solution of a very well known problem of two- 
dimensional electrostatics. Namely, it is easy to see that 0o is the potential 
created by metallic rod of the length n inserted in the uniform electric field. 
We have for the gradient of the field 



d x (j) -id T (f) = — ( 1 



z 2 — n 2 



The action of the optimal configuration 0o is given by a convergent integral 



h 2 



27m 2 7m 2 



Q 2 2Q 2 



Finally, we obtain the leading contribution to the probability of uniformly 
ferromagnetic string P n 



P n ~e~ an , (B.8) 



7T 



a =w (B ' 9) 

We see that a for the string with uniform magnetization is much bigger^] than 
the one for a free string (with only total magnetization fixed) and correspond- 
ing probability decays much faster. 

The solution (B.7) is diverging near the ends of the string. This divergence is 
milder than the one for a free string (c^o ~ r~ 1//2 instead of r _1 ). Therefore, 
we expect bosonization to work even better in this case. Indeed, the estimate 
for r (<9 M O ( r o) ~ ^FL) gives r ~ c 2 n. Therefore, the short distance cutoff is 
defined as 

ro = max jl, c 2 n j . 

We see that for large strings n ^> cr 2 the action coming from the r vicinities 
to the ends of the string is ~ (ttRtq) 2 = c 2 n 2 /Q 2 <C an 2 and the small 
parameter justifying the use of bosonization is c 2 (compare to 1/lnc -1 in the 
case of a free string). Therefore, if one is trying to interpolate (illegally) our 
results to the case of the maximally ferromagnetic (c ~ 1) string one should 
use (B.8,B.9) instead of (38). 

It is also easy to obtain from (B.4) (f(x) = x) the solution O f° r f° r the 
optimal configuration in the limit of high temperature (very long strings). In 
this limit the main contribution to the action comes from 4>o{ z ) ~ 2x/Q when 
< x < n so that S(<f>o) ~ 77? which gives 



4 This is expected because we minimize the action in the smaller class of configu- 
rations. 
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Pn~e" 7n , (B.10) 

Q 1 



7=t£- (B.11) 



It is not a surprise that this result is identical to (40) obtained for the free 
string. When string is very long the main contribution comes from the part 
of the solution which is linearly growing with distance and is uniform in time 
direction. There is also a difference with the case of a free string in that 
the crossover from the Gaussian to the exponential behavior of P n occurs at 
n ~ (3 (not n ~ /51nc _1 as for the free string). The reason is the much milder 
divergence of the solution in the case of a uniform string. 
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